Weak localization of short pulses in disordered waveguides 
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We consider the phenomenon of weak locahzation of a short wave pulse in a quasi- ID disordered 
waveguide. We show that the long-time decay of the average transmission coefficient is not purely 
exponential, in contradiction with predictions of the diffusion theory. The diffusion theory breaks 
down completely for times exceeding the Heisenberg time. We also study the survival probability 
of a quantum particle in a disordered waveguide and compare our results with previous calculations 
using the super-symmetric nonlinear sigma model. 



Since the prediction of Anderson yj that sufficiently 
strong disorder can block propagation of waves and lead 
to localization of wave energy in space, extensive efforts 
have been made to observe this 'Anderson localization' 
for waves of various nature (Schrodinger, acoustic, elec- 
tromagnetic, etc. waves) 0. One particular experiment 
that makes localization evident is sending a short wave 
pulse into a disordered medium and then observing its 
evolution in course of time. If the disorder is weak, the 
wave will propagate out from the source by diffusion (at 
least at distances exceeding the mean free path £) and 
its energy will be eventually distributed over the entire 
space, being negligible at any given point. If, in con- 
trast, the disorder is sufficiently strong, the scattering 
will prevent the wave from going away from the source 
and the energy of the wave will remain localized within 
a volume of linear size ^ around the source. The length 
scale 5 is called the localization length. The experiment 
to 'detect' localization seems then trivial: just wait long 
enough and examine which of the above scenarios is real- 
ized. In reality, however, the situation is complicated by 
the fact that real disordered samples are finite in size and 
therefore, localized or not, the wave leaks out from the 
sample through its boundaries, making the 'final' (i.e. 
corresponding to very long times) state of the system 
identical for both weak and strong disorder |^. This 
implies that in order to distinguish between diffuse and 
localized regimes, one has to analyze the leakage itself, 
and not only the final state of the system. A correct anal- 
ysis of the leakage requires careful treatment of bound- 
ary conditions. Indeed, the leakage happens through the 
boundaries and hence the existence of the latter cannot 
be neglected. 

The purpose of this paper is to present a relatively 
simple theoretical model that captures the main features 
of the diffusion-localization transition correctly account- 
ing for boundary conditions at the surface of disordered 
sample. We then apply the model to describe the 'weak 
localization' phenomenon, a precursor of Anderson lo- 
calization that can be observed in nominally diffusive 
disordered samples. We limit ourselves to the case of 
the so-called 'quasi-lD' experimental geometry that is 
frequently encountered in experiments and is extensively 



studied theoretically. A quasi-lD sample is an open cylin- 
drical tube (waveguide) of length L £, diameter d such 
that X < d < £ {X is the wavelength), and base surface 
A = nd"^ /A. The small diameter of the tube allows one to 
neglect the transverse variation of the average intensity 
of diffuse wave, largely simplifying the derivations. The 
tube has reflecting walls but open ends, and it is flUed 
with a disordered medium (e.g. a mixture of transpar- 
ent and scattering balls in the case of microwaves). A 
short wave pulse is emitted at a position z' inside the 
tube at time t = Q and the intensity of the wave is mea- 
sured at some other position z at time t. In typical ex- 
periments, a wave is incident on the tube from outside 
and the transmitted wave is measured. This corresponds 
to z' ~ £ (since € is a typical distance needed to con- 
vert the incident ballistic wave to the diffuse wave inside 
the sample) and z = L. The following parameters are 
commonly used when light scattering in quasi- ID dis- 
ordered waveguides is considered: number of transverse 
modes N = k'^A/An > 1 (fc 2tt/X) 01, dimension- 
less conductance g = {A/3)N£/L and localization length 
f = {2/3)N£ > £. 

We now consider a short wave pulse incident on the 
waveguide at z = (in this case we set z' — £) or pro- 
duced inside the waveguide (any z' between and L) at 
time t — 0. To find the average transmission coefficient 
T{t) of the tube measured at z = L, we apply the self- 
consistent approach developed in Ref. yl and based on 
the self-consistent theory of localization |j| . It amounts 
to write a diffusion equation for the intensity Green's 
function C{z, z', il): 
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dz dz 



C{z,z',n)^5{z-z') (1) 



supplemented with a self-consistent equation for the 
position- and frequency-dependent diffusion coefficient: 
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D{z,n) Do ^ 
and with the boundary conditions: 



c{z,z,n) 



(2) 



C{z, z', n) ~ zo^^^^^^Ciz, z',n) = Oatz = (3) 
Do dz 
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C{z, z', n) + zo^^^^^^C{z, z', f^) = at z = i (4) 
Do dz 

Here Dq = vi/3 is the 'bare' value of the diffusion co- 
efficient {v is the transport velocity) and zq ^ ^ is the 
so-called 'extrapolation length' that allows one to ac- 
count for reflection of diffuse waves at the sample bound- 
aries. Equation ||2Jl describes the renormalization of 
the diffusion coefficient due to the interference of time- 
reversed trajectories inside the disordered medium. This 
renormalization has the same physical origin as the phe- 
nomenon of coherent backscattering Q. Note that the 
renormalized diffusion coefficient D(z, il) appears not 
only in the diffusion equation ^ but in the boundary 
conditions El as well. 

If 5 ^ oo (or, equivalently, if i/^ 0), waves prop- 
agate by diffusion and the second term on the right- 
hand side of Eq. |5J can be neglected. We then have 
D{z,il) = Dq, and by solving the diffusion equation ^ 
with the boundary conditions © we find 



transmission coefficient: 



Coiz,z',n) = 



L 
Wo 



X [sinh(7z</i) + 7(zo/L) cosh(7z</L)] 

X [sinh(7(f - Zy/L)) + jizo/L) cosh(7(I - z>/L))] 

X [(I + ^\zo/Lf) sinh7 + 27(zo/L) cosh 7] (5) 

where we use the subscript '0' to denote the case of bare 
diffusion, z< — min(z, z'), z> = max(z,z'), and 7 = 

i-iniyDo)^^^. 

We now consider very long waveguides (L » £) and 
neglect the extrapolation length zq ^ £ <^ L in Eq. 
This yields 

L sinh(7z</£)sinh[7(l - z>/L)] 

Coz,z,ri = — — 6 

Dq 7 smh 7 

The average transmission coefficient is then found as 

fom = -Do±Coiz^L,z',n)-'^^^^^ (7) 
az sinn7 
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^ Res \fo{n) ejip{-int),n = -in^TT^Da/L^ 

00 

- ^(-I)"+insin(n7rz7L) exp(-n t/to) (8) 

where to = L'^/t^'^Dq is the time of wave diffusion 
through the waveguide. To evaluate the integral in Eq. 
(0 we used the fact that the integrand To(r2) exp(— ifli) 
has simple poles on the imaginary axis of the complex 
plane at = —in'^Tr'^Do/L'^ and applied the residue the- 
orem. Using Eqs. and (jSJ we can find the steady-state 



Toit)dt = foin = 0) 
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(9) 



(10) 



At long times t > tjj only the first term contributes 
appreciably into the sum of Eq. © and the average 
transmission coefficient decays exponentially with time: 
To{t) ~ {2TTDo/L'^)sm{TTz'/L)exp{-t/tD)- This expo- 
nential decay is a hallmark of diffusion behavior and it is 
often used to test for diffusion behavior of waves in ex- 
periments. The diffusion constant Dq can be extracted 
from the measured average transmission as 

L^d_ 
^dt 

However, recent experiments 0] show that even in the 
diffuse regime {g ^ 1) deviations from the simple expo- 
nential decay of the average transmission coefficient can 
be detected. This phenomenon is called weak localiza- 
tion and it has also an impact on the steady-state {fl = 0) 
transmission coefficient of a disordered waveguide. In the 
latter case, however, the phenomenon is difficult to ob- 
serve because one has to introduce some additional factor 
(e.g. magnetic field) to brake the time-reversal symme- 
try, thus preventing the renormalization of the diffusion 
constant, to be able to compare transmission coefficients 
in the presence and in the absence of the time-reversal 
symmetry. We now show how our theoretical model can 
be used to describe the weak localization phenomenon 
for short wave pulses. 

As far as g remains much larger than unity, we can 
use 1/g as a small parameter and only keep the terms of 
order 1/g (and drop the terms of order 1/g^, 1/g^, ■ ■ ■) 
in Eqs. ^ and (|2Jl. This amounts to write 



D{z,n) 


= Do + -D,{z,n) 
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(11) 


c{z,z\n) 


= Ca{z,z',n) + -Ci{z,z',n) 
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(12) 


f{n) 


- fo{n) + -fl{n) 
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(13) 


m 


= Ta{i) + -Ti{t) 
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(14) 



Substituting these into Eqs. Q and (0) and collecting 
the terms of order 1 / g we find 



Di{z,n) 



.Dt 



Co(z,z,f^) 
Ci(z,z',f7) 



d_ 

dz 



Diiz,n)—Coiz,z\n) 



(15) 



(16) 



We now note that Eq. llt)|) with the delta-function 6{z — 
z') on the right-hand side is also obeyed by Co{z, z', J7). 
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Co{z, z', r2) is therefore the Green's function of Eq. (fTB|l . 
The solution of the latter can then be written as 



Ci{z,z',n)^-4^ I dz"Co{z,z",n) 



' dz" 



Co(z", z", Q)—Co{z", z', n) 



(17) 



The integral in Eq. (|17|l can be evaluated, and for t > t]j 
we obtain 



m = Toit) 



l -9 

Do g 



ao(z') + ai{z')^ + ^2 ( 7^ 



ai{z') + 2a2T— 
to 



(18) 
(19) 



where 



3 / z' 



1 



a2 



and £'(t) is defined by analogy with Eq. (|10|) : 



(20) 
(21) 
(22) 

(23) 



Equations H18|l and (|19() are the main results of this pa- 
per. They lead us to the following important conclusions. 
First, the decay of the average transmission coefficient 
T{t) is not purely exponential and hence, strictly speak- 
ing, the propagation of waves in a disordered waveguide 
does not obey diffusion laws. This is manifested by the 
time dependence of the diffusion coefficient (|19|l which 
has to be time-independent for the purely diffusive prop- 
agation. Second, the diffusion model with a constant, 
time-independent diffusion coefficient can be a good ap- 
proximation only for short times t <^ gto- It is easy 
to show that tu = gto is actually a characteristic time 
scale that is known as the 'Heisenberg time' in the field 
of quantum chaos. The Heisenberg time is the typical 
time needed for a wave (or quantum particle) to visit the 
whole system (the whole waveguide, in our case). Af- 
ter the Heisenberg time, the wave has to pass by those 
parts of the sample that it has already visited before. 
The wave trajectory then crosses itself. Such crossings 
are precisely the phenomenon that renormalizes the dif- 
fusion coefficient and that is accounted for by the sec- 
ond term in Eq. Apparently, our theoretical model 



does not allow us to treat the situations when such cross- 
ings become dominant. We cannot therefore treat times 
exceeding the Heisenberg time. Another way to under- 
stand this limitation of our theoretical approach is to note 
that the Heisenberg time is also the inverse of the typi- 
cal spacing between the modes of a closed sample. For 
times exceeding tn the mode structure of the wave field 
becomes important (i.e. one resolves individual modes). 
This mode structure is not included in our model and 
hence we cannot treat this regime correctly. 

It is worth-noting that ao{z') and ai{z') in Eqs. lfTH|) 
and H19|) depend on the position z' of the source of waves 
and can even change sign when the source position is 
changed. For example, when the wave is incident on the 
waveguide from outside (as in the experiments of Cha- 
banov et al. [^), we set 2' ~ £ and find ao — —3/4 
and ai — — 3/27r^ for L ^ £. On the other hand, if the 
source is placed in the middle of the waveguide, we find 
ao = —7/16 and ai — 3/27r^. In contrast, the value of 
a2 is not sensitive to the source position. This illustrates 
that for long times the waves explore the whole sample 
and have only a limited memory about the position of 
their source. 

The decrease of the diffusion coefficient with time pre- 
dicted by Eq. p9|l has been recently observed in mi- 
crowave experiments by Chabanov et al. (Oj. Out theo- 
retical model allows a reasonably good description of the 
experimental results . It is also worthwhile to note that 
a result similar to Eq. (|18|l has been obtained by Mirlin 
who studied the 'survival probability' P{t) of a quan- 
tum particle in a quasi- ID disordered waveguide using 
the super-symmetric nonlinear sigma model and found 



lnP(.).-(l 



1 t 
TT^g to 



(24) 



for orthogonal symmetry (/3 = 1), corresponding to 
the case of preserved time-reversal symmetry considered 
here. To calculate P{t) in the framework of our self- 
consistent approach we note that the reason for the sur- 
vival probability to differ from 1 is the leakage of wave 
energy through sample boundaries. This leakage is quan- 
tified by the transmission coefficient T(t) and we there- 
fore find 



dP{t) 
dt 



-2f{t) 



(25) 



where the bar denotes averaging over the position of the 
source z' and the numerical factor 2 is the number of 
(identical) open boundaries. For bare diffusion [g — > 00) 
we substitute Eq. into Eq. (|25|l and using the 'initial' 
condition P(0) = 1 find 



-Po(i) = ^exp ( --^ 
TT- V to 



and 



-lnPo(t) = — -ln( — 



(26) 
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£oT t > to. As expected, the decay of Po{t) is purely 
exponential. Note that instead of using the condition 
P(0) = 1 we could arrive at the same result (|26|l by 
requiring P(oo) = 0. 

Consider now finite but large g. Substituting Eq. l(TH|) 
into Eq. (|25|) and performing the integration yields P{t) 
up to an additive constant. The condition P{0) — 1 
cannot be used to determine this constant (as for g — > oo) 
because Eq. (fH^ only holds for t > and does not apply 
for t = 0. We therefore use the condition P{oo) = to 
determine the additive constant and obtain 



Pit) 




t 

ao + «!- h a2 



to 



(27) 



where otQ = Zj-K^ — 3/8, a\ 
This is equivalent to 



= 5/27r2, and aa = I/tt^ 



InP(i) = 



to 



In 



1 t 

2'!T-^g TT^g to 
1/3 3 
9 



(28) 



If we drop the time-independent terms in the square 
brackets [which seems to be done in Eq. (|24|l ] , the resem- 
blance between Eqs. and becomes evident. In- 
deed, the numerical coefficient —l/n'^g in front of {t/tn}'^ 
is the same in both cases. On the other hand, the coeffi- 
cients in front of t/tD differ by 5/2n^g. It is interesting to 
note that the coefficients in front of {t/toY and t/to do 
not depend on the exact condition [P(oo) — 0, P(0) = 1 
or P{t ^0) = Po{t ^ 0)] that we employ to determine 
the additive constant when integrating Eq. (|25|l . 

In conclusion, we have presented a theoretical model 
that enables us to describe the phenomenon of localiza- 



tion for short wave pulses in disordered media, correctly 
accounting for boundary conditions on the sample sur- 
face. The model is based on a diffusion equation with a 
self-consistently renormalized diffusion coefficient. Using 
this model, we calculate the time-dependent transmis- 
sion of waves through a disordered quasi- ID waveguide 
and the survival probability of a wave (or quantum parti- 
cle) in the waveguide in the diffuse regime (dimensionless 
conductance much larger than unity). We show that the 
decay of these quantities with time is not purely expo- 
nential, even for nominally diffusive samples, and slows 
down due to the phenomenon of weak localization. We 
compare our result for the survival probability with that 
obtained previously using a different approach, and show 
that the two results are similar, but not completely iden- 
tical. 
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